A regularized nonnegative canonical polyadic decomposition algorithm with preprocessing for 3D fluorescence spectroscopy

We consider blind source separation in chemical analysis focussing on the 3D fluorescence spectroscopy framework. We present an alternative method to process the Fluorescence Excitation-Emission Matrices (FEEM): first, a preprocessing is applied to eliminate the Raman and Rayleigh scattering peaks that clutter the FEEM. To improve its robustness versus possible improper settings, we suggest to associate the classical Zepp’s method with a morphological image filtering technique. Then, in a second stage, the Canonical Polyadic (CP or Candecomp/Parafac) decomposition of a nonnegative 3-way array has to be computed. In the fluorescence spectroscopy context, the constituent vectors of the loading matrices should be nonnegative (since standing for spectra and concentrations). Thus, we suggest a new NonNegative third order CP decomposition algorithm (NNCP) based on a non linear conjugate gradient optimisation algorithm with regularization terms and periodic restarts. Computer Nadège Thirion-Moreau is the corresponding author. This work has been supported by a “PRES euro-mediterranean” grant, and by the European Research Council under the European Programme FP7/20072013, Grant Agreement no.320594, DECODA project.